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1. Introduction 

O , Since tiie appearance of the pioneering paper of Snyder and Weinberg, many ar- 

' tides have been pubhshed on the problem of meson fields in the presence of strong 

electric and gravitational fields. In 1940 Schiff, Snyder and WeinberJ2l carried out 
one of the earliest investigations of the solution of the Klein-Gordon equation with a 
^ ' strong external potential. They solved the problem of the square well potential and 

discovered that there is a critical point Vcr wh ere t he bound antiparticle mode ap- 
pears to coalesce with the bound particle. PopoJ^Elghowed that this phenomenon is 
a particular effect characteristic of short range potentials and, consequently it should 
not be expected to be observed when one deals with Coulomb interactionJ^EEEl 
The asymptotic limit of the Schiff-Snyder effect, for infinite walls, has been studied 
by Fulling. In 1979, Bawin and LavinJ^ demonstrated that the antiparticle p- wave 
bound state arises for some conditions on the potential parameters, showing in this 
way that, for short range potentials, the presence of the Schiff-Snyder effect strongly 
depends on the angular momentum contribution. 

The presence of strong fields introduces quantum phenomena, such as super- 
criticallity and spontaneous pair production, that cannot be described using per- 
turbative techniques. The discussion of overcritical behavior of bosons requires a 
full understanding of the single particle spectrum, and consequently of the exact 
solutions to the Klein-Gordon equation. Overcritical behavior associated with scalar 
particles in strong short range potentials could be of interest in understanding quan- 
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turn effects such as superradiancJ^^^ or particle production in the vicinity of black 
holes IHEl. 

In the present article, we solve the Klein-Gordon equation for a Woods-Saxon 
and cusp potentials. The interest in computing bound states and spontaneous pair 
creation processes in such potentials lies in the fact that they possess properties 
that could permit us to determine how the shape of the potential affects the pair 
creation mechanism. We show that the antiparticle bound states arise also for the 
Woods-Saxon potential well, which is a smoothed out form of the square well. We 
also show that a cusp potential well supports antiparticle bound states. 

The article is structured as follows: Sec. 2 is devoted to discuss the Klein-Gordon 
equation. In Sec. 3 we solve the Klein-Gordon equation in the presence of a one- 
dimensional Woods-Saxon potential well. We derive the equation governing the 
eigenvalues corresponding to the bound states and compute the bound states. In Sec. 
4 we compute the bound states of the Klein-Gordon equation in the presence of an 
one-dimensional cusp potential well. We also show the dependence of supercritical 
states on the strength and shape of the potential. Finally, in Sec. 5, we briefly 
summarize our results. 



2. The Klein-Gordon equation 

The Klein-Gordon equation minimally coupled to a vector potential A^^ can be 
written as 

//^"(a^ -I- ieAf,){d, + ieA,)^ + m^c^ = 0, (1) 

where the metric ry'^" has signature —2, and we have set h = c = 1. Choosing to 
work with a vector potential of the form, 

eA° = Vir), 1=0, (2) 

Eq. Q can be written as: 
d 



— +iV{r)j (/) - V^</> + m^(/) = 0. (3) 

Eq. Q can be c omple tely separated in spherical coordinates when one works with 
radial potential^^^l Analogously, we have that in Cartesian coordinates, the 
Klein-Gordon equation can be reduced to a second order ordinary differential equa- 
tion when the potential depends only on a single space variable x. In this case Eq. 
© takes the simple form: 

-f [{E V{x)f - rn' - pi] c^{.) = 0, (4) 

where, since the potential V{x) does not depend on time, we have separated vari- 
ables in the form 



0(.T, t) = 4>{x) e:sq){~iEt + ipyU + ipzz), 



(5) 
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Fig. 1. The Woods-Saxon potential well for L 
line) . 



2 with a = 10 (solid line) and a = 3 (dotted 



a.ndp'i ^pI+pI. 

Eq. 10} has the same structure than the one obtained when one solves the Klein- 
Gordon equation for s states, in spherical coordinates, for radial potentialJ^^I 
therefore the results reported in this article are also valid for s-waves in radial 
Woods-Saxon and cusp potentials. In the next section we solve Eq. Q in the pres- 
ence of a Woods-Saxon potential well. 



3. The Woods-Saxon well 

The Woods-Saxon potential is defined a m 

Qi-x) 



V{x) 



-Vo 



e(x) 



(6) 



where Vq is real and positive for a well potential; a > and L > are real and 
positive. Q(x) is the Heaviside step function. The parameter a defines the shape of 
the barrier or well. The form of the Woods-Saxon potential is shown in Fig. ^ In 
order to solve the one-dimensional Klein-Gordon equation 10} , with p±^ = 0, in the 
presence of the Woods-Saxon potential lO we split the solutions into two regions. 
Let us consider bound states solutions for x < 0. In this case we need to solve the 
ordinary differential equation 



dx^ 



E 



Vo 

-_, — a{x+L) 



=0, 



(7) 



where in Eq. |(7} and thereafter we have equated the mass m to unity. Since the 
potential (|S} depends only on one space variable, the contribution of the transverse 
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momentum in Eq. Q appears via the introduction of an additive term that can 
be included in an effective mass term mefj of the form niefj = {m? +p^)^/^. For 
small values of p_L, the contribution of the transverse momentum does not introduce 
qualitative changes to the energy spectrum. 

On making the substitution ~ \ + ^ Eq. becomes 



a^y{\ - y) — 
dy 



y(i - y) 



d4>L{v) 



dy 



{E + V^yY-\ </)L(y)=0. 



(8) 



Setting tphiv) = u'^i^ ^ vT'^iy) and substituting it into Eq. we obtain the 
hypergeometric equation 

y{l - y)h" + [(1 + 2cr) - 2(a + 7 + l)y]h' - (i + d + 7 + A) (i + cr + 7 - A) = 0, (9) 

where the prime denotes a derivative with respect to y and the parameters a, 7, 
and A are 



(10) 



(11) 



a 2a 
The general solution of Eq. (O can be expressed in terms of Gauss hypergeometric 
functions ^1 



h{y) ^ai2Fi[^+-f + a-X,^+-f + a + \,l + 2a;y 



'a2y 



so 



2F1 [ ^ + -f - a - X,^ + J - (7 + XA - 2a;y 



(12) 



^L{y) = aiy^ {1 ~ yy 2F1 Q + 7 + a - A, ^ + 7 + d + A, 1 + 2a; y 

+ a2y-^{l - y)^ 2F1 Q + 7 - a - A, ^ + 7 - (T + A, 1 - 2d; . (13) 
Now we consider the solution for a; > 0. We solve the differential equation 



dx^ 



E 



Mx) = 0. 



After making the substitution z ^ = 1 + e"'^ ^\ Eq. (|14|l can be written as: 



a^z{l — z) — 
dz 



z{l-z) 



d(j)R{z) 



dz 



(14) 



(15) 



Introducing <f)R{z) = 2*^(1 — z) ^g[z) and substituting it into Eq. (fT^ . we obtain 
that 4>r{z) satisfies the hypergeometric equation: 

z(l - z)g"+[{l + 2a) - 2(a- 7+ 1)^]^'- (i + ^- 7 + A) (i + (7-7 - A) 5 = 0, (16) 
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where the prime denotes a derivative with respect to z. The general solution of Eq. 
(PI isEl 

g{z) = feiaFi Q - 7 + cr - A, i - 7 + cr + A, 1 + 2cr; z 

+ b2z'^^2Fi Q-7-fT-A,i-7-(7 + A,l-2(T;z) , (17) 

Mz) = - z)"^ 2i^i Q - 7 + - A, i - 7 + <7 + A, 1 + 2(7; z 

+ 52Z-"(l-z)-T2Fi Q-7-a-A,i-7-a + A,l-2a;z^ . (18) 

As X — > — oo we have that y goes to zero, Analogously, as a; ^ oo, z also goes 
to zero. We choose the regular wave functions 

<l>L{y) - anf{l -yy2Fi(^^+j + (7-X,^+-f + a + \,l + 2a:y 

0fl(z) = 6iz"(l-z)-T2Fi Q-7 + ,T-A,i-7 + 'T + A,l + 2a;z^ . (19) 

In order to find the energy eigenvalues, we impose that the right and left wave 
functions and their first derivatives must be matched at a; = 0. This condition leads 
to 



1 + 2cr 



Mi 



Mi 



+ (i-7+CT-A)(i-7+cr+A) X 



-T + o-+A,2 + 2cr,(l + e 



2^'l(5-7 + o--A,5-7 + <'+A,l + 2o-,(l + s-'»-f')-l) 

+ 2ff(l+e-"-^) = 0, (20) 

which is the eigenvalue condition for the energy E. Explicit solutions of Eq. H20|l . 
giving E in terms of Vq, can be determined numerically. We consider the range 
— 1 < i? < 1 for the values of E. Some aspects of the dependence of the spectrum 
of bound states on the potential strength Vq are shown in Figs. 2 and 3. At some 
value of Vo, a bound antiparticle state appears, it joins with the bound particle 
state, they form a state with zero norm at Vq = Vcr and then both vanish from the 
spectrum. 

The normalization of the wave functions in Eq. (|19|l is given by 



= 2 



dx[E - V{x)](l){x)*(P{x) 



(21) 



The norm of the Klein-Gordon equation vanishes at Vcr, where both possible 
solution and S^") meet. Fig. 2. shows that for 2.0900 < Vq < 2.0908 two 
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states appear, one particle (solid line), and one antiparticle state (dashed line). In 
Fig. 3 the same behavior is observed for 2.3462 < Vq < 2.3463. Particle bound 
states (-E'"*"^) and antiparticle bound states (E^^^) correspond to iV > and < 
respectively. For N = both solutions meet and have the same energy. Antiparticle 
states appear in all the cases considered. For L — 2, we moved the shape parameter 
a from 1 to 18 and, for i = 1, we considered a = 10. Fig. 4 and Fig. 5 show 
the behavior of the turning point (E) versus the potential parameters a and Vq 
respectively. Fig. 4 shows that as the value of a increases, the energy value, for 
which antiparticle states appear, increases. In Fig. 5 we observe that, as the value 
of Vo increases, the energy value for which antiparticle states appear decreases. This 
behavior indicates that square well potentials exhibit antiparticle bound states for 
values of E larger than for smoothed out potentials. 
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Fig. 2. Energy of the lowest bound-state spectrum for L = 2, a = 10. Inset is an enlargement of 
the critical area, showing solid and dotted lines corresponding to positive and negative norm state 
solutions respectively. The critical value for Vq = 2.0908 corresponds to i? = —0.993698. Energy 
is given in units of the rest energy mc^ 



4. The Cusp potential well 

In this section we are interested in studying the Klein-Gordon eq uati on |@J, with 
p± = 0, in the presence of a cusp potential given by the expression^! 

f —Vne^/'^ for T < 

The form of the potential (|22|l is shown Fig. © . The parameter Vq > determines 
the depth of the cusp potential well. From Fig. © one readily notices that the 
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Vo 

Fig. 3. Energy of the lowest bound-state spectrum for L = 2, a = 2. Inset is an enlargement of 
the critical area, showing solid and dotted lines corresponding to positive and negative norm state 
solutions respectively. The critical value for Vq = 2.3463 corresponds to i? = —0.999943. Energy 
is given in units of the rest energy mc^ 
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20 



Fig. 4. Turning point versus a for L = 2. Energy is given in units of the rest energy mc^ 

potential becomes sharper as the shape parameter a becomes smaher. In fact, the 
asymptotic limit of V{x) as a ^ is {-Vo/a)d{x). The potential (|22l) can also be 
regarded as a limit case of the Woods-Saxon potential well as the shape parameter 
a ^ 1 and L = in Eq. ©. As for the Woods-Saxon potential well, we split the 
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Fig. 5. Turning point versus Vq for L = 2. 



solutions into two regions. 
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Fig. 6. The Cusp potential well for a = 2 (solid line) and a = 0.5 (dotted line). 
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Consider the bound states solution for a; < 0. We solve the differential equation 



E + Voe 



X I a 



- 1 



0L(a;) -0. 



On making the change of variables y — 2iaVbe^/'^, Eq. (|23|) becomes 



^ dy \ dy 

-1/2 



(m£; + y/2) +0^ 



Putting (pL ~ y f{y) we obtain the Whittaker differential equation 



f{y)" 



1 iaE 1/4- a2(l - £;2) 



f{y) = 0. 



(23) 



(24) 



(25) 



Ay 

The general solution of Eq. (|25|) can be expressed in terms of Whittaker functions 
M,,^(2/) and W^,,^(y) asHSI 

My) = ciy~'/^M^^iy) + C2y-^'^W,^,{y), (26) 

where ci and C2 are arbitrary constants, and k and are 

K = -iaE, fi = a\/l -E^. (27) 

Analogously, we proceed to obtain the solutions of Eq. Q in the presence of the 
potential H22|l for a; > 0. In this case we have the differential equation 



d^(t)R{x) 
dx^ 



(E + V^e-'/'^y 



cj)R{x) = 0. 



On making the change of variables z — 2iaVoe Eq. (|28|) becomes 



dz V dz 



(iaE + z/2) 



= 0. 



Putting — z ^/^g{z) we obtain the differential equation 



1 iaE 



l/4-a2(l-S2) 



g{z) = 0, 



(28) 

(29) 

(30) 
(31) 



whose general solution is 

0fl(z) = drz-^l''M^^{z) + d2Z-^I^W^^{z), 

where di and ^2 are arbitrary constants. 

Since we are looking for bound states of Eq. with the potential H22|l . we 
choose to work with regular solutions 4>L{y) and 4>r{z) along the x axis: 

(j)L{y) = ciy~^/^M^^^{y), 
Mz) = diz-^/^M,Jz). 

Imposing the condition that the left <t>L{y) and right (t>R{z) scalar wave functions 
must be continuous at a; = as well as their first derivatives, we obtain that the 
energy eigenvalues must satisfy the equation: 

(1 - 2iaVo + 2K)M,^{2taVo) - (1 + 2k + 2fi)M^+i,^{2iaVo) - 0, (33) 



(32) 
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The explicit solutions of Eq. (|33|l . showing the dependence of the energy £' on Vb 
and a can be determined numerically. The dependence of the ground state spectrum 
on the potential strength Vg is shown in Figs. Q and JSJ. We can readily see that 
for a = 0.5, antiparticle states arise for V > 3.6050 and, for a — 1.2, antiparticle 
states arise for V > 3.04386. Here the solutions coming from the lower continuum 
into the region of bound states meet with solutions from the upper continuum. This 
leads to the possibility of spontaneous productions of pairs. As the Woods-Saxon 
potential well, with the cusp potential well, the particle bound states (-E*^"*'-') and 
antiparticle bound states (-E^"-*) correspond to iV > and iV < respectively. 
For = both solutions meet and have the same energy. The normalization of 
the wave functions H32(l is given by Eq. H21|l . Antiparticle bound states arise for all 
values of the shape parameter considered. We varied the shape parameter a from 
0.05 to 2.5. Figs. © and H10|l show the behavior of the turning point (E) versus 
the potential parameters a and Vq respectively. Fig. Q shows that as the value 
of a increases, the energy value, for which antiparticle states appear, decreases. In 
Fig. (|l()|l we observe that, as the value of Vq increases, the energy value for which 
antiparticle states appear increases. 
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Fig. 7. Energy of the lowest bound-state spectrum for a = 0.5. Inset is an enlargement of the 
critical area, showing solid and dotted lines corresponding to positive and negative norm state 
solutions respectively. The critical value for Vq = 3.605 corresponds to E = —0.984766. Energy is 
given in units of the rest energy mc^ 
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Vo 

Fig. 8. Energy of the lowest bound-state speetrum for a = 1.2. Inset is an enlargement of the 
critical area, showing solid and dotted lines corresponding to positive and negative norm state 
solutions respectively. The critical value for Vo = 3.043865 corresponds to E = —0.999996. Energy 
is given in units of the rest energy mc^ 
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Fig. 9. Turning point versus a. Energy is given in units of the rest energy mc^ 

5. Conclusions 

In this article wc have shown that the Woods-Saxon and the cusp potential well are 
able to bind scalar particles. 
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Fig. 10. Turning point versus Vq. Energy is given in units of the rest energy mc^ 

The Woods-Saxon potential, analogous to the square well potential, shows an- 
tiparticle bound states. The turning point, where the norm is zero, depends on the 
potential parameters a and Vq. Therefore the one-dimensional Woods-Saxon po- 
tential exhibits a behavior characteristic of short range potentials 13. It should be 
expected that, like in the Coulomb problem, for slow damping potentials no an- 
tiparticle bound states appear. The results reported in this article suggest that the 
one-dimensional SchifF-Snyder- Weinberg effect El can be extended to the case of 
potentials with non compact support, provided they exhibit, for large values of the 
space parameter, a fast damping asymptotic behavior. 
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